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Abstract 

Micromachined  solid-oxide  fuel  cells  ((jlSOFCs)  are  among  a  class  of  devices  being  investigated  for  portable  power  generation.  Opti¬ 
mization  of  the  performance  and  reliability  of  such  devices  requires  robust,  scale-dependent,  design  methodologies.  In  this  first  analysis,  we 
consider  the  structural  design  of  planar,  electrolyte-supported,  (JtSOFCs  from  the  viewpoints  of  electrochemical  performance,  mechanical 
stability  and  reliability,  and  thermal  behavior.  The  effect  of  electrolyte  thickness  on  fuel  cell  performance  is  evaluated  using  a  simple 
analytical  model.  Design  diagrams  that  account  explicitly  for  thermal  and  intrinsic  residual  stresses  are  presented  to  identify  geometries 
that  are  resistant  to  fracture  and  buckling.  Analysis  of  energy  loss  due  to  in-plane  heat  conduction  highlights  the  importance  of  efficient 
thermal  isolation  in  microscale  fuel  cell  design. 
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1.  Introduction 

There  is  growing  interest  in  the  application  of  microscale 
devices  (also  called  microelectromechanical  systems, 
MEMS)  for  portable  power  generation  [1-11].  It  is  envi¬ 
sioned  that  devices  with  overall  dimensions  on  the  order  of 
millimeters,  and  containing  components  with  critical  dimen¬ 
sions  in  the  range  0.1-10  pun,  will  provide  between  0.5  and 
50  W  of  power  to  enable  the  operation  of  other  microsys¬ 
tems  (such  as  sensors  and  actuators)  and  portable  electronic 
devices  (such  as  radios  and  cell  phones).  The  device  con¬ 
cepts  being  investigated  include  miniaturized  gas-turbine 
engines  [4],  micro  batteries  [5],  proton  exchange  membrane 
fuel  cells  [6-9],  and  solid-oxide  fuel  cells  (SOFCs)  [9-11]. 
The  last  mentioned  is  the  focus  of  this  paper. 

The  principle  of  operation  of  solid-oxide  fuel  cells  is 
well  known  [12,13].  The  active  electrochemical  element  is  a 
three-layer  composite  structure  consisting  of  two  electrodes 
separated  by  a  solid-oxide  electrolyte.  The  essential  fuel 
cell  reaction  is  the  combination  of  hydrogen  and  oxygen 
to  form  water.  The  electrons  participating  in  this  reaction 
flow  from  the  anode  to  the  cathode  through  an  external 
circuit,  generating  direct  electric  current  in  the  process,  and 
the  electrical  circuit  is  completed  by  the  flow  of  oxygen 
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ions  from  the  cathode  to  the  anode  through  the  solid-oxide 
electrolyte.  Such  structures  can  be  realized  in  both  planar 
and  tubular  geometries  [13]. 

The  fabrication  and  performance  of  planar  solid-oxide 
fuel  cells  containing  thick-film  (>10  ptm)  layers,  derived 
from  powder-  and  solution-processing  routes,  has  received 
considerable  attention  [13].  Since  the  ionic  resistance  of 
the  electrolyte  scales  with  the  thickness,  it  is  believed  that 
reducing  the  component  thickness  can  lead  to  improved 
performance.  This  belief  has  motivated  recent  efforts,  us¬ 
ing  a  combination  of  vapor-phase  thin  film  deposition  and 
silicon  micromachining  techniques,  to  manufacture  mi¬ 
croscale  solid-oxide  fuel  cells  (p-SOFCs)  [9-11].  Indeed, 
Jankowski  et  al.  [9]  have  measured  the  performance  of  a 
ptSOFC  containing  a  1.2-pjn  thick  yttria- stabilized  zirconia 
(YSZ)  electrolyte,  0.5-ptm  thick  nickel  anode,  and  0.8-p.m 
thick  silver  anode;  this  composite  membrane  was  supported 
on  a  0.5-mm  thick,  micromachined,  single-crystal  silicon 
platform  (Fig.  1). 

Optimization  of  the  performance  and  reliability  of  mi¬ 
croscale  solid-oxide  fuel  cells  requires  the  development  of 
robust,  scale-dependent,  design  methodologies.  The  ulti¬ 
mate  goal  is  to  identify  combinations  of  materials,  struc¬ 
tural  geometries,  and  processing  routes  that  will  optimize 
the  electrochemical  performance,  minimize  thermal  losses, 
and  ensure  structural  stability  and  reliability.  The  primary 
challenges,  from  the  viewpoint  of  structural  design,  are 
enumerated  below. 
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Fig.  1 .  (a)  Schematic  side-view  of  a  three-layer  solid-oxide  fuel  cell  con¬ 
sisting  of  porous  electrodes  and  dense  electrolyte,  (b)  Critical  dimensions 
and  operating  temperatures.  The  in-plane  dimension,  b ,  is  ~mm,  while 
the  thickness,  h,  is  ~p.m. 

(i)  The  benefits  offered,  and  constraints  imposed,  by  mi¬ 
crofabrication  techniques  are  considerably  different 
from  those  of  conventional  machining  methods.  In  par¬ 
ticular,  micromachining  permits  the  parallel  fabrication 
of  arrays  of  geometrically-complex,  two-dimensional, 
structures,  but  extensions  to  three-dimensions  are 
mostly  limited  to  extrusion-type  processes  [14].  Planar 
designs  of  the  type  shown  in  Fig.  1  are  considerably 
easier  to  realize  compared  to  tubular  structures. 

(ii)  Thin  film  materials  are  invariably  subjected  to  intrin¬ 
sic  residual  stresses  (ayes)  originating  from  the  details 
of  the  growth  processes  [14],  It  is  well  established  that 
intrinsic  stresses  are  strong  functions  of  the  processing 
conditions.  Moreover,  systematic  changes  in  a  few  pro¬ 
cessing  parameters  (argon  gas  pressure  and  bias  voltage 
in  sputter  deposition,  for  instance)  can  lead  to  stress 
changes  of  several  hundred  megapascals  in  metallic 
and  oxide  thin  films.  This  ability  to  tune  the  stress 
by  controlling  the  process  parameters  implies  that  it  is 
both  reasonable  and  essential  to  treat  ny-es  as  a  design 
variable. 

(iii)  Solid-oxide  fuel  cells  operate  at  temperatures  of  sev¬ 
eral  hundred  degrees  centigrade.  This  requires  that  the 
structural  stability  and  reliability  of  a  three-layer  com¬ 
posite  structure  be  guaranteed  over  a  wide  temperature 
range  (say,  20-650  °C).  The  electrode  and  electrolyte 
materials,  and  the  supporting  silicon  structure,  can 
possess  significantly  different  thermomechanical  prop¬ 
erties,  resulting  in  large  thermal  stresses  (oyh)  during 
thermal  excursions.  Failure  due  to  fracture  (in  ten¬ 


sion)  or  buckling  (in  compression)  represents  a  serious 
reliability  concern. 

(iv)  Operation  at  high  temperatures  requires  efficient  ther¬ 
mal  isolation  schemes  to  prevent  heat  loss  from  the 
membrane  device.  The  difficulty  of  this  task  increases 
with  the  degree  of  miniaturization  since,  in  general, 
heat  losses  increase  with  decreasing  device  dimen¬ 
sions.  Furthermore,  it  is  necessary  to  consider  coupled 
thermal  and  structural  design  since  the  structure  that 
minimizes  heat  losses  may  not  be  mechanically  stable. 

1.1.  Scope  of  analysis 

The  preceding  discussion  makes  clear  that  a  natural  start¬ 
ing  point  for  microscale  solid-oxide  fuel  cells  is  the  planar 
design  of  Fig.  1.  In  this  first  analysis,  we  present  simple, 
analytical  models  to  evaluate  the  effect  of  geometry  on 
performance  and  reliability.  Section  2  presents  closed-form 
expressions,  containing  two  free  parameters,  to  estimate  the 
cell  potential  and  power  density  of  a  planar  (xSOFC.  These 
parameters  are  evaluated  using  the  experimental  results  of 
Jankowski  et  al.  [9],  and  the  results  used  to  determine  the 
effect  of  electrolyte  thickness  on  fuel  cell  performance. 

The  electrodes,  which  are  usually  fabricated  using  silver, 
nickel,  or  nickel- YSZ  cermets,  are  porous  (or  contain  mi- 
crofabricated  through-holes)  to  permit  gas  transport  to  the 
electrode/electrolyte  interfaces.  The  yttria- stabilized  zirco- 
nia  electrolyte  layer  is  dense  and  relatively  thick,  and  serves 
as  the  load-bearing  structure.  Therefore,  as  a  first  approx¬ 
imation,  the  structure  can  be  modeled  as  a  single  -layer 
membrane  with  material  properties  corresponding  to  those 
of  YSZ.  Section  3  analyzes  the  mechanical  stability  and  in¬ 
tegrity  of  this  membrane  and  presents  a  design  diagram  to 
identify  geometries  that  are  resistant  to  both  buckling  and 
fracture.  Heat  losses  due  to  in-plane  conduction,  and  the 
need  for  coupled  design  considerations,  are  discussed  in 
Section  4. 


2.  Electrochemical  performance 

Consider  the  fuel  cell  structure  shown  in  Fig.  1.  The  tem¬ 
perature  of  the  membrane  is  assumed  to  be  uniform  (i.e. 
To  =  7],).  The  primary  performance  parameters  of  the  fuel 
cell  are  the  cell  potential  (E)  and  power  density  (P)  (i.e.  the 
power  produced  per  unit  membrane  area),  which  are  related 
by 

PQ)  =  i  ■  E(i )  (1) 

where  i  is  the  current  density.  In  general,  the  cell  potential 
can  be  expressed  as  [15] 

E  —  E  o  —  TOhmic  —  Activation  —  ^Concentration-  (2) 

The  open  circuit  voltage,  Eo,  is  dependent  solely  on  the 
thermodynamics  of  the  fuel  cell  reaction  and  has  a  value, 
given  by  the  well  known  Nemst  equation  [13],  of  ~1.0V 
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for  typical  operating  conditions  under  consideration.  How¬ 
ever,  when  current  is  drawn  from  the  cell,  the  cell  potential 
decreases  due  to  several  sources  of  polarization,  rj,  as  dis¬ 
cussed  below. 

One  source  of  polarization  is  due  to  the  ohmic  resistivity 
of  the  solid-oxide  electrolyte  to  the  flow  of  oxygen  ions  and 
is  given  by 

*7  Ohmic  =  i~  (3) 

°e 


where  he  and  rxe  are  the  thickness  and  ionic  conductivity 
of  the  electrolyte,  respectively.  The  ionic  conductivity  of 
YSZ  ranges  from  0.3  to  3.0  S/m  over  the  temperature  range 
600-1000  °C,  respectively. 

Both  activation  and  concentration  polarizations  are  related 
to  the  kinetics  of  fuel  cell  reactions.  The  first  is  due  to  the 
activated  processes  involved  in  charge  transfer  at  the  elec¬ 
trode/electrolyte  interfaces  and  can  be  expressed  as  [15] 


Here  R  is  the  gas  constant  (8.3 14  J/K/mol),  F  the  Fara¬ 
day  constant  (96485  C/mol),  n  the  number  of  electrons 
exchanged  per  oxygen  ion  (hence,  n  =  2),  and  To  the 
operating  temperature  of  the  fuel  cell  (assumed  to  be  uni¬ 
form  over  the  membrane).  The  exchange  current  densities 
i\  and  T  are  typically  treated  as  free  parameters  and  ob¬ 
tained  from  fits  to  measurements.  Each  of  the  two  terms  in 
Eq.  (4)  corresponds  to  the  limiting  cases  of  the  more  general 
Volmer-Butler  analysis  of  the  kinetics  of  activated  processes 
[18]. 

The  concentration  polarization  is  due  to  diffusion-limited 
processes  (such  as  transport  of  fuel  to  the  active  sites 
on  the  membrane)  in  the  operation  of  the  fuel  cell  [15]. 
Analytical  expressions  for  the  concentration  polariza¬ 
tion  of  conventional  fuel  cells  are  usually  complicated 
functions  of  microstructural  parameters  of  the  electrodes 
and  electrolytes.  In  this  first  analysis,  it  is  assumed  that 
^Concentration  =  0.  This  is  reasonable  since,  the  smaller  the 
dimensions  of  the  device,  the  smaller  the  relevant  diffusion 
lengths,  and  the  less  important  any  concentration-induced 
polarizations. 

With  this  assumption,  and  combining  Eqs.  (l)-(4),  we 
obtain 


P  =  iE  =  i\E0-  - 

L  < 


-a)]' 


These  expressions  are  independent  of  the  shape  of  the  mem¬ 
brane  since  the  temperature  is  assumed  to  be  uniform. 

The  exchange  current  densities  are  evaluated  by  fitting 
(5)  and  (6)  to  the  measurements  of  Jankowski  et  al.  [9]. 


Fig.  2.  Electrochemical  performance  of  microscale  solid-oxide  fuel  cell. 
The  data  correspond  to  the  measurements  of  Jankowski  et  al.  [9].  The 
curves  correspond  to  an  analytical  model  that  treats  the  exchange  current 
densities  as  fitting  parameters. 


Fig.  2  shows  that  values  of  200  and  450  A/m2  result  in  good 
agreement  with  the  data  obtained  at  600  °C.  These  values 
are  comparable  to  the  exchange  current  densities  reported 
in  thick-film  SOFCs  at  similar  temperatures  [19].  Values  for 
exchange  current  densities  at  490  and  550  °C  were  obtained 
by  similar  analysis  (Table  1).  This  table  shows  that  the  av¬ 
erage  power  density  in  ixSOFCs  is  ~103  W/m2  at  the  tem¬ 
peratures  indicated.  Hence,  a  power  output  of  1 W  requires 
a  circular  membrane  of  radius  ~2  cm. 

Eqs.  (5)  and  (6)  can  be  used  to  evaluate  the  effect  of 
electrolyte  thickness  on  performance.  Fig.  3  indicates  that, 
even  for  membranes  as  thick  as  2p,m,  the  activation  po¬ 
larization  dominates  over  ohmic  losses  at  the  temperatures 
considered,  and  a  further  reduction  in  thickness  is  not  ex¬ 
pected  to  result  in  improved  power  density.  Therefore,  the 
choice  of  membrane  thickness  can  be  guided  primarily  by 
mechanical  and  thermal  considerations  as  explained  in  the 
following  sections. 


Fig.  3.  Comparison  of  activation  and  ohmic  polarizations  for  a  microscale 
solid-oxide  fuel  cell. 
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Table  1 

Measured  and  fitted  p.SOFC  parameters 


Temperature,  To 

763  K  (490  °C) 

823  K  (550  °C) 

873  K  (600  °C) 

Ionic  conductivity  of  YSZ,  o-e  (S/m)  [17] 

0.1 

0.16 

0.3 

Peak  power  density,  P  (W/m2)  [9] 

360 

730 

1450 

RTqIF 

0.065 

0.071 

0.075 

Exchange  current  density,  i\  (A/m2) 

75 

150 

450 

Exchange  current  density,  12  (A/nr) 

50 

100 

200 

3.  Structural  reliability  and  stability 

The  composite  membrane  shown  in  Fig.  1  is  invari¬ 
ably  subjected  to  intrinsic  (trres)  and  thermal  (trth)  residual 
stresses.  It  is  well  known  that  layered  structures  are  suscep¬ 
tible  to  stress-induced  failure  by  several  modes  including 
fracture,  buckling,  interlaminar  delamination,  and  spalbng 
[20].  Examining  all  these  modes  in  a  composite  membrane 
complicates  the  design  process  significantly.  As  mentioned 
previously,  however,  we  consider  here  a  class  of  structures 
in  which  two  layers  are  porous  and  the  third  (electrolyte) 
is  relatively  thick,  dense,  brittle,  and  serves  as  the  primary 
load-bearing  structure.  Hence,  for  first-order  design,  the  de¬ 
sign  of  this  single  membrane,  assumed  to  be  circular  with 
radius  b  and  thickness  h,  against  fracture  and  buckling  is 
investigated.  The  material  is  assumed  to  be  isotropic  and 
homogeneous,  and  has  a  Young’s  modulus  E,  Poisson’s 
ratio  v,  and  coefficient  of  thermal  expansion  a.  (The  ther¬ 
mal  expansion  coefficient  of  the  supporting  structure  is  as; 
Table  2.)  As  before,  it  is  assumed  that  To  =  7],. 

Since  the  membrane  is  supported  by  a  relatively  large  and 
stiff  silicon  structure,  it  is  reasonable  to  assume  clamped 
boundary  conditions.  The  validity  of  this  assumption  was 
verified  using  finite-element  computations.  Similarly,  since 
the  gaseous  streams  on  either  side  of  the  membrane  are 
at  approximately  equal  pressures,  we  assume  that  there  is 
no  pressure  differential  across  the  membrane.  The  thermal 
stress,  which  can  be  readily  obtained  by  integrating  the  ther¬ 
mal  strain  over  the  diameter  of  the  membrane  and  subse¬ 
quently  employing  Hooke’s  law,  is  given  by  [21] 

^(a  -  as)(7b  -  Trt)  =  A“  AT’  (7) 

where  Trt  is  room  temperature.  A  positive  AT,  which  cor¬ 
responds  to  expansion  of  an  unclamped  membrane,  induces 


Table  2 

Nominal  material  properties  at  300  K  [16] 


Yttria-stabilized 
zirconia  (YSZ) 

Silicon 

Young’s  modulus,  E 

200  GPa 

160  GPa 

Poisson  ratio,  v 

0.2 

0.2 

Coefficient  of  thermal 

10  x  10“6  K”1 

3  x  10-6K-' 

expansion,  a 

Density,  p 

6000  kg/m3 

2320  kg/m3 

compressive  thermal  stresses  in  the  clamped  membrane 
(crth  <  0),  as  expected. 

Estimation  of  intrinsic  residual  stresses  is  non-trivial  and 
usually  requires  experimentation,  especially  for  the  complex 
oxide  materials  of  interest  in  fuel  cells.  Previous  measure¬ 
ments  suggest  that  it  is  possible  to  tune  the  intrinsic  stress 
in  YSZ  thin  films  over  a  large  range  by  suitable  choice  of 
processing  technique  (sputtering  and  electron-beam  evapo¬ 
ration,  for  instance)  and  processing  parameters  [11,22-24], 
It  is  assumed  here  that  crres  is  independent  of  temperature, 
although  this  is  strictly  true  only  if  the  film  microstructure 
remains  invariant  during  thermal  excursion. 

3.1.  Stress-induced  failures 

The  total  stress  (a)  in  the  membrane  is  the  sum  of  the 
thermal  and  intrinsic  components  and  can  be  expressed  as 
E 

(4  =  °res  +  CTth  =  ffres  “  - A  a  AT.  (8) 

1  -  V 

If  the  stress  in  the  membrane  is  sufficiently  compressive, 
then  failure  will  occur  via  buckling.  For  a  circular  plate  that 
is  clamped  at  the  edges,  the  critical  stress  at  which  buckling 
occurs  is  given  by  [25,26], 


Thus,  large  area,  thin  membranes  (i.e.  small  ( h/b )),  which 
are  desirable  from  the  viewpoint  of  electrochemical  perfor¬ 
mance,  are  seen  to  be  particularly  susceptible  to  buckling. 

The  primary  mode  of  failure  in  tension  is  brittle  fracture 
due  to  the  nucleation  and  growth  of  cracks  [27].  The  simplest 
criterion  assumes  that  fracture  occurs  when 

ff  >  fff  (10) 

where  erf  is  the  fracture  strength  of  the  structure.  However, 
the  fracture  strength  of  brittle  materials  (such  as  yttria- 
stabilized  zirconia)  is  a  distributed  quantity  and  probabilis¬ 
tic  methods  are  required  for  reliable  design.  In  the  widely 
used  Weibull  method,  such  design  is  accomplished  by  first 
fitting  experimentally  measured  fracture  strengths  of  nomi¬ 
nally  identical  specimens  (of  volume  Vo)  to  a  two-parameter 
Weibull  function  [27].  These  parameters  are  the  reference 
strength,  op,  which  is  the  stress  at  which  63%  of  the  speci¬ 
mens  fail,  and  a  Weibull  modulus,  m,  which  is  a  measure  of 
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the  degree  of  scatter  in  the  fracture  strength.  The  probability 
of  failure  of  a  structure  of  volume  V  can  be  expressed  as 

(id 

where  eri,  <T2,  and  a 3  are  the  principal  stresses  in  the  struc¬ 
ture.  The  stress  state  in  the  membranes  can  be  approximated 
to  be  equibiaxial  such  that  03  =  rr2  =  o'  and  <73  =  0.  For 
the  purpose  of  design,  it  is  convenient  to  specify  a  tolera¬ 
ble  probability  of  failure,  and  compute  the  maximum  stress 
that  can  be  imposed  on  the  structure.  Using  the  geometry 
defined  in  Fig.  1,  and  using  Eq.  (11),  the  maximum  stress 
can  be  expressed  in  terms  of  the  Weibull  parameters  as 

<i2> 

The  Weibull  parameters  for  vapor-deposited  YSZ  thin  films 
are  currently  unknown. 

3.2.  Design  diagram 

Eqs.  (8),  (9),  and  (12)  can  be  used  to  construct  a  design 
diagram  (Fig.  4),  with  axes  corresponding  to  the  intrinsic 
residual  stress  and  the  imposed  temperature  change.  This 
map  shows  the  effects  of  these  parameters  on  the  total  stress 
in  the  membrane.  The  stresses  are  compressive  above,  and 


tensile  below,  the  reference  line  corresponding  to  a  =  0. 
The  stress  limits  for  buckling  and  fracture  are  shown  on  the 
map.  A  further  limitation  is  placed  by  the  practical  consid¬ 
eration  that  AT  >  0  (i.e.  the  membrane  is  typically  never 
cooled  below  room  temperature  either  during  fabrication  or 
operation)  and  that  the  membrane  be  stable  at  AT  =  0.  Us¬ 
ing  these  limits,  a  safe  design  space  can  be  highlighted  on 
the  map.  Clearly,  the  maximum  permissible  temperature  ex¬ 
cursion  can  be  increased  by  increasing  the  membrane  aspect 
ratio  (h/b).  The  same  goal  can  also  be  achieved  by  decreas¬ 
ing  Aa,  although  the  limited  choice  of  electrochemically 
active  materials,  and  substrates  for  microfabrication,  is  ex¬ 
pected  to  restrict  this  option. 

With  this  approximation,  and  for  membranes  with  neg¬ 
ligible  intrinsic  residual  stress,  Fig.  5  shows  a  plot  of  the 
maximum  permissible  temperature  excursion  as  functions  of 
membrane  geometry.  For  a  1-pm  thick  membrane  operating 
at  600  °C,  this  graph  suggests  a  maximum  in-plane  radius  of 
only  17  pm.  This  in  itself  does  not  represent  a  serious  lim¬ 
itation  since  large  arrays  of  small,  circular,  membranes  can 
be  processed  in  parallel  to  obtain  the  requisite  large  (~cm2) 
active  areas. 

Alternately,  the  mechanical  stability  of  the  membrane 
can  be  increased  either  by  choosing  other  shapes  or  by 
using  stiffening  structures.  Fig.  6  shows  the  computed  crit¬ 
ical  buckling  stresses  for  clamped,  equal-area,  membranes 
of  different  shapes  and  indicates  that  rectangles  have  mod¬ 
estly  higher  resistance  to  buckling  compared  to  squares 
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In-plane  Membrane  Radius,  b  (m) 

Fig.  5.  Maximum  permissible  temperature  change  for  yttria-stabilized 
zirconia  membranes  as  functions  of  radius  and  thickness.  The  intrinsic 

and  circles.  A  more  promising  strategy  is  to  pattern  silicon 
nitride  stiffeners  as  supports  for  the  zirconia  film,  as  illus¬ 
trated  in  Fig.  7.  Preliminary  finite-element  calculations  in¬ 
dicate  that  the  effect  of  such  stiffeners  is  to  impose  clamped 
boundary  conditions  on  each  sub-unit  of  the  electrolyte 
membrane,  thereby  increasing  the  resistance  to  buckling. 
Detailed  numerical  analysis  of  the  mechanical  behavior  of 
stiffened  membranes  is  an  area  of  our  current  research. 


4.  Thermal  performance 

Thermal  losses  may  occur  during  the  operation  of  fuel 
cells  through  convection  to  the  working  fluids,  radiation  to 
the  surroundings,  and  conduction  in  the  structure.  Conduc¬ 
tive  losses  are  considered  here  since  they  directly  influence 
the  structural  design  of  the  device.  In  contrast  to  the  as¬ 
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Fig.  6.  Normalized  critical  buckling  stress  for  equal-area  rectangular 
membranes  as  functions  of  the  aspect  ratio  of  side  lengths.  The  circles 
indicate  the  corresponding  stress  values  for  circular  membranes  of  the 

sumption  of  uniform  operating  temperature  of  the  previous 
sections,  we  now  consider  the  membrane  to  be  subjected  to 
an  in-plane  gradient  such  that  the  temperature  is  To  at  the 
center,  and  7j,  =  7rt  at  the  boundaries.  Practically,  such  a 
profile  may  be  maintained  by  distributing  thin-film  heaters 
on  the  membrane  and  cooling  the  substrate. 

Under  the  assumption  of  uniform  heating,  the  temperature 
distribution  in  the  membrane  as  a  function  of  radial  position, 
r,  is  parabolic  [28], 

T(r)  =  7o  -  (7o  -  7b)  Q2  .  (13) 

The  power  per  unit  membrane  area,  Q,  required  to  maintain 
this  temperature  profile  is 

e  =  *(p)<T°“Thi  =  K?)A7''  (14) 


Fig.  7.  Schematic  of  patterned  stiffeners  that  may  be  used  to  prevent  buckling  of  thin-film  fuel  cell  membranes.  Stiffeners  offer  the  possibility  to  increase 
the  buckling  load  without  significantly  increasing  the  thermal  conduction  losses  from  the  membrane. 
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curves.  Heat  loss  increases  with  decreasing  thickness  because  thinner 
membranes  require  a  smaller  diameter  to  resist  buckling. 

where  k  is  the  thermal  conductivity  of  the  membrane.  The 
thermal  analysis  can  now  be  coupled  with  the  preceding 
structural  considerations  by  combining  Eqs.  (7)-(9).  Assum¬ 
ing  zero  intrinsic  stress,  we  obtain 

AT2 

Q  =  3.28(1 -v) - k  Act  (15) 

h 

Eq.  (15)  represents  the  thermal  power  density  required  to 
maintain  a  membrane,  which  is  required  to  be  resistant 
to  buckling ,  at  temperatures  Tq  and  7rt  at  the  center  and 
periphery,  respectively.  This  analysis  can  conveniently  be 
presented  in  graphical  form,  with  separate  vertical  axes 
corresponding  to  the  critical  membrane  radius,  b,  and  Q, 
plotted  against  the  temperature  change,  AT.  Contrary  to  in¬ 
tuition,  Fig.  8  indicates  that  the  heat  losses  increase  as  the 
membrane  is  made  thinner.  This  can  be  understood  on  the 
basis  that,  although  making  the  membrane  thinner  reduces 
conduction  heat  losses,  it  also  requires  a  corresponding 
decrease  in  the  membrane  radius  to  maintain  structural 
stability.  In  other  words,  this  coupled  analysis  captures  a 
scaling  behavior  that  cannot  be  obtained  by  examining  the 
thermal  performance  or  the  structural  behavior  in  isolation. 

For  a  given  area,  a  circular  shape  has  the  smallest  perime¬ 
ter  and  hence  the  lowest  conduction  heat  loss.  Even  for  this 
optimal  shape,  it  is  striking  to  note  that  heat  loss  densities 
are  considerably  larger  than  the  generated  power  density  for 
temperature  gradients  as  small  as  tens  of  degree  centigrade. 
Developing  efficient  thermal  isolation  schemes  for  p,SOFCs 
is  therefore  critical  for  ensuring  efficient  fuel  cell  operation. 


5.  Summary 

Microscale  solid-oxide  fuel  cells  are  among  a  class  of 
devices  being  investigated  for  portable  power  generation. 


As  a  step  toward  developing  systematic,  scale-dependent, 
design  methodologies  for  such  devices,  we  presented  ana¬ 
lytical  models  for  the  electrochemical  performance,  struc¬ 
tural  reliability  and  thermal  losses  for  a  planar,  circular, 
electrolyte-supported,  device.  The  principal  assumptions  and 
results  are  enumerated  below. 

(i)  Expressions  for  the  cell  potential  and  power  density 
were  presented  based  on  the  assumption  that  concen¬ 
tration  polarizations  are  negligible  in  microscale  fuel 
cells.  Values  for  the  exchange  current  densities  were 
obtained  by  fitting  these  expressions  to  previously  mea¬ 
sured  fuel  cell  performance.  These  results  were  then 
used  to  evaluate  the  relative  importance  of  ohmic  and 
activation  polarizations.  Our  results  indicate  that  fuel 
cell  performance  is  insensitive  to  electrolyte  thickness 
in  the  micrometer  range  at  the  temperatures  considered. 

(ii)  Layered  composite  thin-film  structures  are  invariably 
subjected  to  intrinsic  and  thermal  residual  stresses.  Fail¬ 
ure  due  to  fracture  and  buckling  are  hence  serious  re¬ 
liability  concerns.  Design  diagrams  were  presented  to 
identify  geometries  that  are  stable  and  reliable. 

(iii)  The  thermal  performance  of  planar  microscale  fuel 
cells  was  analyzed  assuming  in-plane  heat  conduction 
to  be  the  dominant  loss  mechanism.  This  analysis, 
coupled  with  structural  considerations,  enables  the 
identification  of  geometries  that  simultaneously  min¬ 
imize  thermal  losses  and  ensure  sufficient  structural 
stability.  However,  over  the  range  of  parameters  con¬ 
sidered,  thermal  losses  are  significantly  higher  than  the 
electrochemical  power  produced.  Developing  schemes 
for  thermal  isolation  of  microscale  fuel  cell  devices  is 
therefore  an  important  area  of  research  in  this  field. 
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